Bistable non-volatile elastic membrane memcapacitor exhibiting chaotic behavior 
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We suggest a realization of a bistable non- volatile memory capacitor (memcapacitor). Its design 
utilizes a strained elastic membrane as a plate of a parallel-plate capacitor. The applied stress 
generates low and high capacitance configurations of the system. We demonstrate that a voltage 
pulse of an appropriate amplitude can be used to reliably switch the memcapacitor into the de- 
sired capacitance state. Moreover, charged-voltage and capacitance-voltage curves of such a system 
demonstrate hysteresis and transition into a chaotic regime in a certain range of ac voltage ampli- 
tudes and frequencies. Membrane memcapacitor connected to a voltage source comprises a single 
element nonautonomous chaotic circuit. 



I. INTRODUCTION 

Currently, there is a strong interest in resistive, capac- 
itive and inductive elements with memory^. The mem- 
ory feature extends functionality of such elements and 
leads to novel circuit application^!^. The main recent 
progress in this area is in the field of memory resistive 
(memristive) systems^ H that can be potentially used in 
both digitaP and analojJMMED applications. 

This article concerns less studied memcapacitive sys- 
tems. By definition 1 , a voltage-controlled memcapacitive 
system is given by the equations 

q(t) = C(x,V C ,t)V C (t) (1) 

x = f(x,V c ,t) (2) 

where q(t) is the charge on the capacitor at time t, Vc(t) 
is the applied voltage, C is the memcapacitance, a; is a 
set of n state variables describing the internal state of the 
system, and / is a continuous rt-dimcnsional vector func- 
tion. It is important that the memcapacitance C depends 
on the state of the system and can vary in time. Several 
systems showing memcapacitive behavior has been iden- 
tified i ncludi ng vanadium dioxide metamaterials^, ionic 
systems^E^ and superlattice memcapacitorsSS. 

In this paper we suggest a different memcapacitor re- 
alization. Our main idea is to replace a plate of a 
parallel-plate capacitor by a strained elastic membrane 
as we demonstrate schematically in Fig. [I] The ap- 
plied stress bends the membrane up or down, allowing 
for two equilibrium positions. When the membrane is 
in the position closer to the bottom plate, the capac- 
itance of the device is higher (for simplicity, this con- 
figuration is called "1"). When the membrane is bent 
up, the system has a lower capacitance denoted by "0". 
Both states are perfectly stables, providing a better non- 
volatile information storage capability than the known 
designs of memcapacitors^"^!. We show below that the 
information can be written into the memcapacitor state 
by applying a single voltage pulse. An interesting feature 
of such a device is a chaotic behavior regime achievable 
within a certain range of parameters. 



II. MODEL 

Let us consider an air gap capacitor composed by a 
fixed plate and a flexible strained membrane (bottom 
and top plates in Fig. [lj. For the sake of simplicity, 
we describe the membrane by a single variable z (the 
effective displacement of the membrane from its middle 
(non strained) position) and calculate the capacitance us- 
ing a parallel plates capacitor model. The effect of the 
stress can be described by a double well potential of the 
form U(z) oc (z 2 — z 2 ) , where ±zo are the equilibrium 
positions of the membrane. In addition, there is an at- 
tractive electrostatic interaction between the charges on 
the fixed plate and membrane. The electric force acting 
on the membrane is given by the product of the charge 
on it, q, times the electric field produced by the charge 
on the fixed plate, —qj (2eo 5), where eo is the electric 
constant and S is the area of the fixed plate. 

Introducing a dimensionless variable y — z/d, where d 
is the separation between the bottom plate and middle 
position of the membrane, we write the relation between 
the charge q(t) and voltage V(t) and the membrane's 
equation of motion as 

q(t) = i^ta v ® = civVWt), (3) 




FIG. 1: Schematic of a bistable non- volatile elastic membrane 
memcapacitor connected to a voltage source V(t). The top 
plate of a regular parallel-plate capacitor is replaced by a 
flexible strained membrane. Because of two equilibrium posi- 
tions of the membrane (represented by a double- well potential 
sketched on the right hand side of the figure) , stable high and 
low capacitance configurations are possible in such a system. 
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FIG. 2: Dynamics of membrane memcapacitor subjected to a single voltage pulse. It is shown in (a) ((b)) that a /3o = 2.1 
(/3o = 2.8) pulse drives the system into the final high (low) capacitance state independently of the initial state of the membrane. 
The final state (at the indicated time) is shown as a function of the pulse amplitude fto = [2tt/ (cjo d)] \J Co / (2 m) Vo in (c). 
The red (blue double-dot-dash) curve represents the final position of the membrane when the initial state is "0" ("1"). Below 
fio — 1.77 the final state depends on the initial one. Above Po = 1.77 both curves overlap and the final state is independent on 
the initial state of membrane. The parameters used in these simulations are yo = 0.2 and F = 0.7. 
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where Co = e Q S/d, yo — Zo/d, Y = 2ir"f/u> , 7 is the 
damping constant, uio is the natur al angular frequency 
of the system, /3(t) = [2tt / (ui d)] ^/C / (2 m) V(t), m is 
the mass of the membrane and the time derivatives in 
Eq. Q are taken with respect to the dimensionless time 
t = two/ (2tt). Eqs. ©and Q have the general form 
given by Eqs. ([I]) and ([2]) and, therefore, the membrane 
memcapacitor is a second-order voltage-controlled mem- 
capacitive system. 



III. INITIALIZATION OF MEMCAPACITOR 
STATE AND HYSTERESIS LOOPS 

A single voltage pulse of an appropriate duration and 
amplitude can be used for the purpose of the memca- 
pacitor state initialization. Specifically, we demonstrate 
below that depending on the pulse amplitude, the mem- 
capacitor can be controllably set into the state " 0" or " 1" . 
The initialization of the memcapacitor state can not rely 
on the electrostatic interaction alone because this inter- 
action in the capacitor is always attractive. Therefore, in 
order to obtain the state "0", a higher amplitude pulse 
should be applied such that the restoring elastic force, 
when the voltage pulse ends, allows overcoming the po- 
tential barrier shown in Fig. [I] 

Figs. [2^a) and [2^b) show simulation results of the 
memcapacitor dynamics obtained as a numerical solution 
of Eq. Q for two different values of the applied pulse 
amplitude. When a lower amplitude /3q = 2.1 pulse is ap- 
plied to the system for a long enough time, the membrane 



first finds an equilibrium position at this value of volt- 
age. This position is independent on the initial state of 
the membrane in the switching regime (/? > 1.77 in our 
simulations, see Figs. [2]ja) and[2^c) for details). Then, 
when the voltage drops to zero, the membrane oscillates 
around —yo and eventually the equilibrium position " 1" 
is reached. 

When a higher amplitude pulse (flo = 2.8) is applied 
(see Fig. [2]Jb)), the membrane is attracted down stronger 
before being released. As a result, when the pulse ends, 
the membrane goes up and stays in the low capacitance 
configuration " 0" . With a further increase of the pulse 
amplitude, the final state of the membrane can be again 
"1" and then "0", etc. Such periodic situation is clearly 
seen in Fig. [2jc). Note that in the switching regime 
(j3o > 1.77) the final state of the membrane does not 
depend on its initial state. 

When the membrane memcapacitor is subjected to ac 
voltage, its dynamics can be periodic or chaotic depend- 
ing on the frequency and amplitude of the ac voltage. 
In this section we consider only the periodic case. The 
chaotic behavior is discussed in Sec. IIVI A numerical so- 
lution of Eqs. ( 3p l with (3(t) — /3 sin(u;i) is plotted in 
Fig. [3] for three different values of the applied voltage 
frequency. The hysteresis curves demonstrate typical be- 
havior of memory elemental non-linear dependencies at 
lower frequencies, pinched hysteresis loops at intermedi- 
ate frequencies and linear behavior at higher frequencies. 

The hysteresis observed in Fig. [3] at w = 0.5cj is ba- 
sically related to the well-known phase shift between pe- 
riodic applied force and position of a driven mechanical 
system. Basically, this effect originates from the mem- 
brane's inertia as well as from the fact that a finite time 
is required to dissipate instantaneous energy stored, e.g., 
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FIG. 3: q—V and C—V (top curves in the top inset) curves for 
a membrane memcapacitor calculated when a sinusoidal volt- 
age is applied. The bottom curve in the top inset represents 
y(fi) at oj — 0.5cjo- Parameters used in this simulation are 
the same as in Fig. [5] and /3 = /3o sin(w£) with /3o = 2.8. The 
graphs demonstrate a periodic steady-state regime beyond an 
initial equilibration interval (r > 30). Limits to linear (or- 
ange) and non-linear (blue) behaviors are obtained at high 
and low frequencies. Hysteresis loops are seen at intermedi- 
ate frequencies. The bottom inset is the energy dissipated by 
the system, given by U — J V dq, calculated at u = 0. Swo- 
rn the elastic degree of freedom of the membrane. Since 
the electrostatic interaction between capacitor's plates is 
always attractive, the C—V curve is symmetric with re- 
spect to = and the membrane experiences a driving 
force of the doubled frequency of ac voltage. The pinched 
behavior additionally proofs that the system is a mem- 
ory capacitive device. The bottom inset in Fig. [3] shows 
the energy dissipated by the device for 10 cycles in the 
steady regime. Because of the damping term in Eq. Q, 
the membrane memcapacitor is a dissipative system. 



IV. CHAOTIC BEHAVIOR 

Most experts would agree defining chaos as the ape- 
riodic, long-term behavior of a bounded, deterministic 
system that exhibits sensitive dependence on initial con- 
ditions and system parameters^. In this section, we 
demonstrate that the membrane memcapacitor described 
by equation ^ can exhibit chaotic behavior. In partic- 
ular, chaotic regime is observed in numerical simulations 
for a specific range of amplitudes and frequencies of the 
applied sinusoidal voltage V = V sin(wi). 

Specifically, in order to explore chaotic properties of 
the system, let us investigate how the state of the mem- 
capacitor (at a long selected moment of time) depends 
on one parameter of the system, such as the ac-voltage 
amplitude Vq. For this purpose, we ran multiple simula- 
tion of the system dynamics at different values of f3 and 
monitor the value of capacitance C at a certain moment 



of time. Fig. R|a) shows the value of capacitance evalu- 
ated at r = ICp as a function of the amplitude /3o- This 
plot demonstrates that at small (/3 < 2.25) and high 
(Po > 2.58) voltages, the system reaches practically the 
same value of capacitance (at a given time) if the voltage 
amplitude is slightly changed. However, at intermediate 
voltages (2.25 < f3 < 2.58) predictability in the value of 
the capacitance is lost, which is an indicator of chaotic 
behavior. 

The Largest Lyapunov Exponent (LLE) A is a quanti- 
tative measure of chaos^. It is defined for a long time 
behavior of a single trajectory and is sensitive to the ini- 
tial conditions of the system. If LLE is positive (A > 0) 
the trajectory is chaotic and A is a measure of the av- 
erage rate at which predictability is lost. On the other 
hand, if A < the evolution is not ch aotic. We use a 
standard approach to calculate the LLE 17 * 18 ! Fig. |4^b) 
shows the result of our calculation of LLE as a function of 
the applied voltage amplitude /3rj. The chaotic behavior 
(positive values of A) is clearly observed in the interval 
2.25 < fa < 2.58. 

The system trajectories in the phase space {y vs. y) at 
low, high and intermediate voltage amplitudes are shown 
in Fig. |4jc)-(e) correspondingly. If the voltage amplitude 
is such that /3 is below approximately 2.25 (Fig. [4jc)) , 
the membrane moves periodically either in the low or 
high capacitance region. If the voltage amplitude is high 
(/3 is above 2.58), a different kind of periodic steady-state 
solution is observed (see Fig. [4jd)) . In this case, the volt- 
age is high enough for the membrane to pass through the 
barrier between the potential wells generating oscillations 
of the membrane between the high and low capacitance 
regions. At both low and high voltage amplitudes, the 
system dynamics is not chaotic. Finally, at intermediate 
voltage amplitudes (2.25 < /3q < 2.58), the system does 
not reach any steady solution. Figs. [4^e-l) and |4je-2) 
show the system trajectories up to two different simula- 
tion times (r = 100 and r = 1000). It is evident that 
the system is ergodic, what is another important mani- 
festation of chaos. In fact, it is not quite surprising that 
Eq. Q produces chaotic trajectories. Generally, Eq. Q 
can be seen as a modification of the Duffing oscillator 
equation, whose chaotic properties are well-known^. 



V. CONCLUSIONS 

The membrane memcapacitor is a modification of a 
parallel-plate capacitor, in which one of the plates is re- 
placed by a strained membrane having two equilibrium 
positions. This feature defines two states ("0" and "1") 
that can be used in memory applications. The system 
is identified as a second-order voltage-controlled memca- 
pacitive system. Pinched hysteresis loops for both capac- 
itance and charge as a function of voltage were obtained. 
In addition, a chaotic behavior was demonstrated in a 
certain range of applied voltage amplitudes. 

Several system parameters were used in our numeri- 
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FIG. 4: Chaotic behavior of the membrane memcapacitor subjected to ac-voltage (/? = fio sin(wt)). The value of capacitance 
(a) at r = 10 4 and the Largest Lyapunov exponent (b) are plotted as a function of the amplitude /3o. The remaining plots 
(c)-(e) demonstrate trajectories in the phase space (y vs. y) of the system at three different values of the amplitude /3o- Plots 
(c) and (d) show a periodic steady-state behavior, and (e) represents the ergodic behavior at /3o = 2.4. Calculation parameters 
are the same as in Fig. [2] and cj — 0.5 u . 



cal simulations. From engineering standpoint, the most 
important parameter is z$ (denning the equilibrium po- 
sitions of the membrane). The choice of zq should pro- 
vide well defined states "0" and "1". Moreover, a con- 
tact between the plates should be avoided at any time. 
Our calculations indicate that both conditions are satis- 
fied at z = 0.2<i Smaller (larger) values of z (at fixed 
values of all other parameters) will decrease (increase) 
the local potential maxima at z = and thus facilitate 
(complicate) the mechanical switching between low- and 
high-capacitance states. Moreover, a variation of zq will 
modify the range of parameters in which the system ex- 
hibits the chaotical behavior. Basically, it is expected 
that at smaller (larger) values of z chaos will be observed 
at smaller (large) applied voltage amplitudes. Similarly, 



the window of non-chaotic behavior will be modified, al- 
though the main frequency features of hysteresis loops 
will remain unchanged. The chaotic behavior will disap- 
peared in the limit zq — > . 

It is important to note that our calculations show that 
chaos is possible in an ac-driven single- element electronic 
circuit. Previously suggested chaotic circuits-EHS! i n _ 
volving those with memristors^, involve several circuit 
elements. Moreover, the membrane memcapacitor is a 
passive device in contrast with the active Chua's diode^ 
and active memristors^l used in existing chaotic circuits. 
Finally, we are not aware about any other capacitive de- 
vices exhibiting chaos. The membrane memcapacitor can 
possibly be fabricated using standard MEMS or NEMS 
(nano-electro-mechanical system) fabrication techniques. 
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